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41 Abstract: Conditions required to support buoyant convective plumes are investigated for 

42 explosive volcanic eruptions from circular and linear vents on Earth, Venus, and Mars. Vent 

43 geometry (linear versus circular) plays a significant role in the ability of an explosive eruption to 

44 sustain a buoyant plume. On Earth, linear and circular vent eruptions are both capable of driving 

45 buoyant plumes to equivalent maximum rise heights, however, linear vent plumes are more 

46 sensitive to vent size. For analogous mass eruption rates, linear vent plumes surpass circular vent 

47 plumes in entrainment efficiency approximately when L a > 3 r 0 owing to the larger entrainment 

48 area relative to the control volume. Relative to circular vents, linear vents on Venus favor 

49 column collapse and the formation of pyroclastic flows because the range of conditions required 

50 to establish and sustain buoyancy is narrow. When buoyancy can be sustained, however, 

51 maximum plume heights exceed those from circular vents. For current atmospheric conditions 

52 on Mars, linear vent eruptions are capable of injecting volcanic material slightly higher than 

53 analogous circular vent eruptions. Flowever, both geometries are more likely to produce 

54 pyroclastic fountains, as opposed to convective plumes, owing to the low density atmosphere. 

55 Due to the atmospheric density profile and water content on Earth, explosive eruptions enjoy 

56 favorable conditions for producing sustained buoyant columns, while pyroclastic flows would be 

57 relatively more prevalent on Venus and Mars. These results have implications for the injection 

58 and dispersal of particulates into the planetary atmosphere and the ability to interpret the 

59 geologic record of planetary volcanism. 

60 

61 1 . Introduction 

62 The mechanisms by which ash and volatiles are transported by explosive volcanism on Venus 

63 and Mars are important issues that relate to present day observations of the surface and the 
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64 


evolution of the climate, atmospheric chemistry and subsurface conditions. Advances in 

65 understanding the range of magma compositions on Mars [ e.g ., Rieder et al., 1997; Bandfield et 

66 al., 2000] continue to expand notions of the prevalence of extensive explosive volcanism. 

67 Detection of elevated SO 2 in the upper atmosphere of Venus [Esposito, 1 984; Moore et al, 1992; 

68 Robinson et al., 1995; Glaze, 1999] and diffuse deposits near the summits of some volcanoes 

69 [ Campbell and Rogers, 1 994] also seem consistent with explosive volcanism. 

70 There is abundant geologic evidence on both Venus and Mars for large volume volcanic 

71 eruptions from long linear vent systems (e.g., Figure 1). Extensive flood volcanism and very long 

72 duration fissure eruptions have been conjectured to have played an important role in the 

73 formation of many units and features on Mars (e.g., Greeley and Spndis [1981]; Keszthelyi et al. 

74 [1999]; Greeley et al. [2000]; Keszthelyi et al. [2006]). It is becoming increasingly evident that 

75 this form of volcanism is also important for Venus. Roughly 20% of the large flow fields on 

76 Venus originate from fractures (e.g., Figure lb) [Stof'an and Smrekar, 2005], and clearly fall into 

77 the flood lava category, with areal extents ranging from 39,000 to 744,000 km 2 [Magee and 

78 Head, 2001], and thickness estimates ranging from <100 m to 1 km [Roberts et al., 1992; 

79 Lancaster et al., 1995], Terrestrial flood basalt eruptions, such as the Roza eruption of the 

80 Columbia River Basalt Group and those associated with the Deccan province, have been 

81 responsible for significant volatile degassing with associated impacts on climate [Thordarson 

82 and Self, 1996; Self et al., 2006]. However, studies that deal with the volcanic contribution of 

83 volatiles into planetary atmospheres [Greeley, 1987; Wilson and Mouginis-Mark, 2001] have 

84 given only limited attention to the ability of plumes associated with explosive phases of large 

85 effusive eruptions to transport and disperse volatiles and ash. 
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86 Buoyant plumes rising above basaltic fissure eruptions on Earth are well documented [e.g., 

87 Thordarson and Self, 1993; Walker et al., 1 984; Lockwood et al., 1 984; Thordarson and Larson, 

88 2007]. Terrestrial experience [ Thorarinson , 1969; Carracedo et al., 1992; Thordarson and Self, 

89 1993 ] indicates that fissure eruptions typically comprise multiple effusive and explosive 

90 episodes. These episodes commonly involve localized activity along various sub-sections of the 

9 1 fissure as the eruption progresses. Active fissure lengths for historical eruptions range from a few 

92 hundred meters (e.g, 1961 eruption of Askja, Iceland [Thorarinson and Sigvaldason, 1962], to a 

93 few kilometers (e.g., 1886 eruption of Tarawera, New Zealand [ Walker et al., 1984], and 1783 

94 eruption of Laki [ Thordarson and Self, 1993]). The surface expression of fissure vents on Venus 

95 and Mars are often tens of km in length. The fissure segment shown in Figure la is just part of 

96 the overall fissure. Concentric fissures, like those at Alba Patera, are an order of magnitude 

97 larger. In addition, volume eruption rates may have been orders of magnitude larger than 

98 terrestrial analogs [ Wilson and Mouginis-Mark, 2001; Baloga et al., 2003; Baloga and Glaze, 

99 2008], Although volcanic fissures visible on the surface of Venus or Mars were not likely active 

100 along their entire length at any one time, large conjectured volumetric eruption rates may point to 

101 active fissure lengths significantly greater than the documented terrestrial analogs. 

102 Over the last 30 years, there have been numerous models for explosive volcanic eruption 

103 columns based on the ideas for buoyant, convecting plumes presented by Morton et al. [1956] 

104 (e.g., Wilson [1976]; Wilson et al. [1978]; Woods [1988]; Stothers [1989]; Thornhill [1993]; 

105 Woods [1993b] Glaze et al. [1997]; Glaze [1999]; Llort and Weitz [2001]). However, with the 

106 exception of Stothers [1989] and Woods [1993a], the effects of linear geometry (including 

107 multiple aligned point sources) on plume dynamics have largely been ignored. In general, 

108 volcanic plume models have focused on eruptions from localized circular vents under the 
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109 assumption that linear vent plumes that rise much higher than the length of active fissure can be 

1 10 approximated with a point source. The original basis for this assumption was the field work by 

1 1 1 Walker et al. [1984] on the 1886 fissure eruption of Tarawera, the largest historical eruption in 

1 12 New Zealand, where most of the erupted material originated from a series of closely space vents 

113 along a 7 km fissure. During its most explosive phase, a Plinian plume rose to an estimated 

1 14 height of 28 km above a 4 km active segment of the fissure. The distal airfall deposit from this 

115 Plinian plume exhibits a symmetry consistent with (i.e., indistinguishable from) fallout from a 

116 plume originating from a circular vent. Walker et al. [1984] infer from the surface distribution of 

117 distal ash that, because the rise height was so much greater than the active fissure length, the 

118 plume material “forgot” that it originated from an extended source. However, this inference 

119 derives from the final surface deposit, and as such, reflects more on the processes that influence 

120 settling rather than vertical transport. It ignores the dynamic transport characteristics that might 

121 be associated with different vent geometries and the ability of a plume to transport material up to 

1 22 a diffusive regime. 

123 It is not at all unexpected that some component of the fallout deposit from a large linear vent 

124 eruption should be more symmetric in nature. First, for narrow linear plumes, most of the 

125 entrainment occurs along the linear extent of the plume, and thus expansion perpendicular to the 

126 fissure axis is dominant. Expansion parallel to the fissure axis should be small in comparison. 

127 Thus, by the time a plume reaches heights such as 28 km, the expansion of the plume may well 

128 be more circular in nature. Second, atmospheric diffusion and turbulence during ash fallout 

129 [Suzuki, 1983; Armienti et al., 1988; Glaze and Self, 1991] produce a randomizing effect on the 

130 final deposition of ash on the surface, smoothing and masking (i.e., “forgetting”) the original 

131 linear nature of the plume rise. However, Walker et al. [1984] also noted that the thickness of the 
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132 proximal airfall deposit did not vary much with distance along the 4 km active segment of the 

133 fissure. From this observation, they concluded that the entire 4 km fissure segment was likely 

134 active concurrently. This observation directly supports the inference that the linear vent 

135 geometry is important to the dynamics of the plume transporting the material found in the 

136 proximal deposit. Thus, dynamics of plume rise, particularly in the first few kilometers, may be 

1 37 significantly influenced by the vent geometry. 

138 Key issues for planetary explosive volcanism include the heights attainable by convecting 

139 plumes, the ability to transport and disperse particulates and volatiles, and the subsurface 

140 conditions within the feeding conduit required to sustain a column. All the models based on 

141 extensions of the Morton et al. [1956] system of equations indicate that atmospheric entrainment 

142 in the first few kilometers of the eruption column play a dominant role in the dynamics of 

143 buoyant plumes. Entrainment, in turn, depends upon the surface area and upward velocity of the 

144 plume in the lower reach. 

145 In this paper, quantitative constraints are developed for the dimensions and character of 

146 buoyant columns for a range of active linear and circular vent systems and eruption conditions 

147 on the Earth, Venus, and Mars. The influence of atmospheric profile, gravity, geometric and 

148 scale differences among these planets are evaluated. This is accomplished by adapting previous 

149 models for buoyant plumes arising from central vents [ Glaze and Baloga, 1996; Glaze et al., 

150 1997; Glaze 1999; Glaze and Baloga, 2002] to linear vent geometries. The models used in this 

1 5 1 work are analogous to the more elementary single-component linear source model of Stothers et 

152 al. [1986] and Stothers [1989] with the adaptation to both solid and vapor plume constituents. 

153 The results obtained in this paper provide significant implications for the relative roles of 
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154 widespread airfall and pyroclastic flow mechanisms in forming the geologic record from 

1 55 explosive volcanism on the terrestrial planets. 

156 2. Convective Plume Model 

1 57 The approach to modeling linear volcanic vents described here builds on the original work by 

158 Stothers [1989] but takes advantage of substantial improvements that have been made in 

159 volcanic plume modeling over the last 20 years. Changes from the Stothers [1989] approach 

160 include (1) explicit inclusion of the conservation of thermal energy, as proposed by Woods 

161 [1988], and implemented by Glaze and Baloga [1996], Glaze et al. [1997], Glaze [1999] and 

162 Glaze and Baloga [2002], (2) consideration of multiple plume components, particularly solids, 

163 and (3) more detailed information on Venus and Mars atmospheric compositions. 

164 In general, the complete system of equations describing buoyant plume rise that must be 

165 solved requires at least a half dozen differential equations and another half dozen equations 

166 describing parameters within the plume and ambient atmosphere. Table 1 shows the system of 

167 differential equations from Glaze et al. [1997] where the control volume and entraining area are 

168 given in a generic form in the first column (Notation is defined in Table 2). For the cylindrically 

169 axisymmetric system in Glaze et al. [1997], shown in Figure 2a, the control volume is defined as 

170 V - m^dz, and the area through which ambient is entrained is A e = 2 w dz. For the linear vent 

171 system shown in Figure 2b, V = 2bLdz and A e ~ 2 Ldz (where it is assumed that L » b). The 

172 approach is then to solve the sets of governing equations in Table 1 for both central and linear 

173 vent geometries. The results described below have assumed that L remains relatively constant as 

174 a function of z. Thus L does not appear in any of the linear vent conservation equations in Table 

175 1. 

176 Boundary conditions at the vent for all variables ( r , b, u, 0, and n) are denoted by a 
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177 subscripted “o”. The vent is defined as the top of the feeding conduit, and the point at which the 

178 erupted material is free to entrain ambient atmosphere and expand. Note that the dimensions at 

179 the top of the conduit can be considerably smaller than the visible expression of the vent at the 

180 surface for both circular and linear vents. Values for the entrainment constant, a, have been 

181 estimated for buoyant plumes originating from both circular and linear vents using empirical 

182 [e.g., Batchelor, 1954; Morton et al., 1956], theoretical [Agrawal and Prasad, 2003], and 

183 simulation [Suzuki and Koyaguchi, 2010] approaches. Values for circular vent plumes range 

184 from 0.05 - 0.15, and from 0.0975 - 0.16 for linear vent plumes. For the purpose of comparing 

185 the effects of geometry only, this study assumes a value of a= 0.09 for all cases. 

1 86 In addition to the conservation equations defined by (3) - (8), the system is completed by the 


1 87 following definitions for volume fraction, 

188 </> d +<t> v +</> s = 1 (9) 

1 89 the ideal gas law, 

190 P = (pJ d R d + pJ v R v )0 (10) 

191 bulk density, 

192 P b =pA+pA+Ps</> s (11) 

193 and the bulk specific heat, 

194 C B =[pJ d C a + pJ v C v + p s <f> s C s ]/p B (12) 


195 Following Glaze and Baloga [1996], Glaze et al. [1997], Glaze [1999] and Glaze and Baloga 

196 [2002], the systems of ordinary differential equations in Table 1 are solved simultaneously using 

197 a Runge-Kutta technique. To avoid numerical artifacts in the boundary region near the plume 

198 top, the maximum plume height is defined as the height at which the upward velocity drops to 
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199 below 10 m/s. Plume rise beyond this point is minimal (relative to the overall plume height), 

200 while the velocity drops quickly and the corresponding plume radius (or half width for linear 

201 plumes) increases exponentially. 

202 

203 3. Results for Earth’s Atmosphere 

204 Much has been published on the sensitivities of the central vent model given in Table 1 [e.g., 

205 Wilson, 1976; Wilson et al., 1978; Woods, 1988; Glaze and Baloga, 1996; Glaze et al., 1997], 

206 Thus, here, the focus is on the primary differences that result from the linear vent geometry 

207 relative to the central vent. The numerical code used to solve the linear vent system is identical to 

208 the central vent model, where the only differences are those resulting from the geometry of the 

209 vents. 

210 Figure 3 illustrates a basic comparison of predicted plume heights on Earth for a range of 

211 vent sizes. As noted above, the length of fissure vent that is active does not appear in the 

212 conservation equations in Table 1. Thus, a simple comparison can be made of plume height as a 

213 function of either b 0 (for linear vents) or r 0 (for circular vents). For the results shown in Figure 3, 

214 all of the boundary values at the vent are kept the same, except for the size of the vent opening 

215 (radius, r 0 , for circular vents, or half-width, b a , for linear vents). The vent conditions used for 

216 both geometries are velocity, u 0 = 300 m s' 1 , temperature, 0 o = 1000 K, and the mass fraction of 

217 water vapor, n 0 = 0.03, where the mass fraction of vapor in the plume is defined as n- p v <j) v ! ' ps- 

218 For this example, the plumes are assumed to rise from sea level into the US Standard atmosphere 

219 (e.g., Wallace and Hobbs [1977]). 

220 For this set of boundary conditions both the linear and circular vent models produce a 

221 maximum plume height of around 35 km. Figure 3 shows maximum plume heights for circular 


Glaze et al. : Buoyant Volcanic Plumes from Linear Vents 


Page 9 



222 


and linear vents, plotted for a broad range of vent sizes. From this plot, it can be seen that 

223 buoyant, convecting plumes originating from circular vents can be maintained with substantial 

224 maximum heights over a wide range of vent sizes, from a radius of a few meters to around 900 

225 meters. For the boundary conditions used in this analysis, plumes from circular vents collapse, or 

226 are otherwise unable to maintain a steady column, for vent sizes greater than about 900 m. 

227 However, the linear plumes are much more sensitive to the vent size, and can maintain a 

228 convective plume over a much more narrow range of vent sizes. For the same boundary 

229 conditions (temperature, velocity, gas mass fraction), plumes from linear vents only reach 

230 significant heights for values of b 0 of a few meters to -200 meters. For this set of boundary 

23 1 conditions, the peak plume height occurs when b 0 ~ 75 m. 

232 Because one of the main drivers of plume height is the mass flux at the vent [Woods, 1988; 

233 Glaze and Baloga, 1996], a better comparison between linear and circular vents can be made by 

234 examining predicted plume heights for equivalent mass fluxes. The mass flux at the vent is 

235 defined as the density of the erupted material multiplied by the cross sectional area of the vent 

236 and the exit velocity (= pb 0 UoA 0 ). Keeping the bulk plume density and exit velocities the same 

237 for both linear and circular vents, the key parameter for comparison becomes the cross sectional 

238 vent area, where A a = 7tr a for the circular vent, and A a = 2 b o L 0 for the linear vent. 

239 Column 1 of Table 3 indicates a range of cross sectional vent areas. Column 2 of Table 3 

240 contains the corresponding radii for circular vents. For linear vents, the half-width, b 0 , 

241 corresponding to the cross sectional area is dependent upon the choice for active fissure length, 

242 L 0 . Noting that the original definition of a linear vent included the assumption that L 0 » b 0 , the 

243 table only considers values for b 0 that are less than an order of magnitude smaller than L 0 . Active 

244 fissure lengths up to 10 km are considered. 
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245 Figure 4 illustrates the plume heights corresponding to the vent sizes given in Table 3. From 

246 Figure 4, it can be seen that as L 0 increases, linear plumes become more capable of establishing a 

247 convective regime over a broad range of b 0 , similar to the circular vents. This is primarily 

248 because as L 0 increases, the entrainment area of the linear plumes increases relative to the control 

249 volume. The ability of a plume to become buoyant is driven by whether or not sufficient air can 

250 be entrained (and warmed) to reduce the bulk plume density before upward momentum is 

251 exhausted. From the right hand side of (3) the mass of dry air entrained by the linear plume is 

252 roughly the same as the cylindrical case when uL « nur. Thus, the linear plumes surpass circular 

253 vents in entrainment efficiency approximately when L 0 > 3 r 0 . Note, however, that because L does 

254 not appear in the governing equations, the maximum plume height in Figure 4 always occurs 

255 around b 0 ~ 75 m. 

256 

257 4. Applications to Venus and Mars 

258 Analysis of the linear vent model under ambient conditions found on Venus and Mars are 

259 now considered. Because of the extreme differences in the atmospheric densities, boundary 

260 condition values used in these analyses are slightly different than those used for Earth. 

261 For comparison with conditions considered by Glaze [1999], the variables held constant 

262 under Venus conditions are the initial velocity, u 0 = 270 m s' 1 , initial temperature, 0 o = 1400 K, 

263 and the mass fraction of water vapor, n 0 - 0.05. Further, the scenarios discussed here have used 

264 the high latitude atmospheric temperature and pressure profile [ Seijf et al., 1985], as explosive 

265 eruptions under these conditions were shown by Glaze [1999] to be the most likely to rise to 

266 significant altitudes. The primary difference between the boundary conditions used here and 
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267 those used by Glaze [1999] is that here it has been assumed that the vent is at 0 km above the 

268 mean planetary radius (as opposed to 1 0 km ampr). 

269 Figure 5 shows the plume heights reached for circular and linear vents for a range of radii 

270 and half widths, respectively. There are several points of interest to note. First, consistent with 

271 work by Thornhill et al [1993] and Glaze [1999], the range of conditions over which buoyant 

272 plumes can be maintained on Venus is very narrow (note the range on the x axis is only 300 m, 

273 compared to > 800 m on Earth). Also, the shape of the curve for circular vents is substantially 

274 different from the terrestrial case, with an almost linear relationship between vent radius and 

275 maximum plume height. The primary explanation for this difference in shape is the much higher 

276 pressure ambient atmosphere. 

277 The maximum plume height achievable for the given boundary conditions is ~52 km 

278 corresponding to a circular vent radius of 247 m. Like the terrestrial case, there is a very abrupt 

279 break at the point where a convective plume can no longer be maintained. For linear vents, the 

280 range of conditions that can maintain a convective plume is likewise very narrow, with the 

281 maximum plume height of ~69 km occurring when b Q = 16 m. Unlike the terrestrial case, linear 

282 vents on Venus appear capable of driving a plume to somewhat higher maximum altitudes, for 

283 all other things remaining equal, albeit under very limited conditions. 

284 As for Earth, a better comparison between the circular and linear vents is made by looking at 

285 the plume heights for each when the mass flux at the vent (and thus the cross sectional vent area) 

286 is the same. Table 4 indicates circular vent radii and linear vent cross widths for a range of cross 

287 sectional vent areas. Again, the linear vent half widths are restricted to those that are less than an 

288 order of magnitude smaller than the fissure length, consistent with the assumption that b 0 « L a . 
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289 Figure 6 shows the predicted plume heights for the vent sizes given in Table 4. It can be seen 

290 that the maximum plume height for the linear vent case occurs consistently for b a = 16 m, but the 

291 peak is shifted to the right by increasing the vent size (i.e., L). Unlike the terrestrial example, 

292 however, the linear vent is more effective than its circular counterpart at producing a convective 

293 plume that is capable of reaching significant heights for the same mass flux. 

294 For Mars, the variables held constant are initial velocity, u 0 = 230 m s' 1 , initial temperature, 

295 6 0 = 1000 K, and mass fraction of water vapor, n 0 = 0.01. It should be noted that Glaze and 

296 Baloga [2002] placed limits on the applicability of buoyancy models such as the one described 

297 here. That study indicated that such models are not valid above heights greater than 10-20 km 

298 on Mars under current atmospheric conditions. For those boundary conditions that produce 

299 plumes higher than about 1 0 km, the alternative approach of Wilson and Head [2007] is used to 

300 model the subsequent inertial transport of particles along ballistic trajectories. Using this 

301 approach, the upward velocity and plume dimensions at 10 km are used as starting conditions for 

302 momentum driven rise. Plume heights found using the convective model are compared to the 

303 heights found using the combined buoyancy-ballistic approach for the circular and linear vents in 

304 Figure 7. 

305 Using the combined buoyancy-ballistic approach, a comparison of the linear and circular 

306 vents is made by examining plume heights for fixed mass fluxes (Table 5, Figure 8), as was done 

307 for Earth and Venus. The ballistic constraint above 10 km results in no plumes that rise above 15 

308 km above the vent. The plumes from circular vents decrease in height as the mass flux (i.e., vent 

309 area) increases, implying column collapse conditions. However, in contrast, the linear vent 

310 plumes continue to show an increase in maximum plume height as mass flux increases, 

31 1 indicating that the upward velocity of the bulk plume material at 10 km altitude (where transport 
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312 transitions from convective to ballistic) is also increasing. The implication is, therefore, that 

313 these linear vent plumes are not necessarily collapsing, and that linear vent plumes develop a 

314 pyroclastic fountain capable of injecting a substantial amount of ash and volatiles into the 

3 1 5 atmosphere at the plume top. 

316 The preference for the formation of pyroclastic flows and near vent deposition on Mars, as 

317 opposed to widely dispersed airfall deposits, is consistent with some prior suggestions in the 

318 literature. Most notably, interpretations of early volcanic activity (e.g., Mouginis-Mark et al. 

319 [1988]; Carr [2007]) have suggested numerous and extensive pyroclastic flows on the flanks of 

320 Alba Patera prior to a transition of late stage effusive volcanism. Similar considerations may 

321 apply to incised units thought to be either airfall or pyroclastic flow deposits at Hecates Tholus 

322 and Hadriaca and Tyrrhena paterae [. Mouginis-Mark et al., 1982; Greeley and Crown, 1990; 

323 Gregg and Farley, 2006]. 

324 

325 5. Conclusions 

326 When buoyancy can be sustained, plumes from both linear and circular volcanic vents have 

327 very similar maximum heights. However, the conditions required to sustain buoyant plumes are 

328 highly restricted for linear vents. The range of conditions for sustained buoyancy is very narrow 

329 for Venus and relatively more favorable for the Earth. The dependence of maximum plume 

330 height on circular vent diameter is dominated by the nature of the ambient atmospheric profiles 

331 for density, pressure, and water content. On Venus, the extreme ambient atmospheric pressures 

332 (100 times Earth) result in an almost linear relationship between the vent radius and plume 

333 height. 
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334 When buoyancy is sustained, linear vents appear to be equally capable of injecting volcanic 

335 ash and volatiles into the atmosphere. For analogous mass flux rates at the vent, the maximum 

336 heights to which linear and cylindrical plumes can rise are essentially equivalent. In cases where 

337 the length of active linear vent is less than approximately 3 r 0 , the entrainment area of linear 

338 plumes is significantly less than the cylindrical (circular vent) equivalent. Thus the range of vent 

339 widths that can sustain a buoyant plume is narrow and these plumes are more likely to collapse. 

340 Only the Earth, with its thick wet atmosphere, favors explosive eruptions that can maintain 

341 convective plumes reaching 10s of km in altitude. Conversely, the range of plausible conditions 

342 that can produce a buoyant convective plume on Venus is very restricted. The probability of an 

343 eruption occurring with just the right boundary conditions is relatively small. Thus, it seems that 

344 pyroclastic flows are much more likely in the nonterrestrial environments. On Mars, the results 

345 presented here favor interpretations of pyroclastic flows and near vent deposition at major 

346 volcanic centers, including Alba, Hadriaca, and Tyrrhena paterae and Hecates tholus, where the 

347 discriminating morphologic diagnostics remain ambiguous. 

348 
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494 Figure Captions 

495 Figure 1. Explosive phases from fissure vents, such as those shown here for (a) Mars (1 8.1° N 

496 245.3° E; THEMIS visible image V05484014) [Mouginis-Mark and Christensen , 2005], 

497 and (b) Venus (5.5° S, 206° E, Magellan SAR image) may have contributed significant 

498 volatiles to the atmosphere. 

499 

500 Figure 2. Schematic views of (a) circular and (b) linear vent geometries. The control volume for 

501 each plume is indicated by the gray shaded volume. The entrainment area is the outer, 

502 exposed surface area of the control volume. The volcanic vent in both cases is defined as 

503 the top of the conduit. Boundary conditions at the vent are denoted by variables with “o” 

504 subscript. Note that for linear vents, L 0 is » b 0 , and that L does not change substantially 

505 as a function of z. 

506 

507 Figure 3. Maximum predicted plume heights as a function of either vent radius (for circular 

508 vents), r a , or half width (for linear vents), b Q . Models assume terrestrial ambient 

509 atmospheric conditions with temperature and pressured defined by the US Standard 

510 Atmosphere, u 0 = 300 m s’ 1 , 0 o = 1000 K, and n 0 = 0.03. 

511 

512 Figure 4. Maximum predicted plume heights as a function of vent area. Model boundary 

513 conditions are the same as those in Figure 3. The linear vent results are shown for 

514 multiple choices of active fissure length. 

515 
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Figure 5. Maximum predicted plume heights on Venus as a function of either vent radius (for 
circular vents), r a , or half width (for linear vents), b a . Models assume northern Venus 
latitude ambient atmospheric conditions, u Q = 270 m s' 1 , 0 o = 1 400 K, and n 0 = 0.05. 

Figure 6. Maximum predicted plume heights on Venus as a function of vent area. Model 

boundary conditions are the same as those in Figure 5. The linear vent results are shown 
for multiple choices of active fissure length. 

Figure 7. Maximum predicted plume heights on Mars as a function of either vent radius (for 
circular vents), r 0 , or half width (for linear vents), b 0 . Models assume the current low 
pressure ambient atmospheric conditions, u 0 = 230 m s' 1 , 9 0 = 1 000 K, and n 0 — 0.01 . (a) 
Predicted plume heights for circular and linear vents when convective model is allowed 
to run to its full extent, regardless of inconsistencies and non-physical effects (e.g., 
expansion rates that exceed the speed of sound), (b) Predicted heights for circular and 
linear vents when convective model is terminated at 10 km. Upward velocity at 10 km is 
then used to estimate the maximum height for a purely ballistic vertical trajectory. 

Figure 8. Maximum predicted plume heights on Mars as a function of vent area. Model 

boundary conditions are the same as those in Figure 7. The convective model has been 
terminated at 10 km (as in Figure 7b) and subsequent rise determined from a ballistic 
vertical trajectory. The linear vent results are shown for multiple choices of active fissure 
length. 
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539 Table 1 . System of generic 1 st order partial differential equations for convective rise, after Glaze 

540 et al. [ 1 997]. The control volume and entrainment area can be defined to accommodate either a 

541 circular or linear source. Notation is defined in the Table 2. 


Description 

Generic Conservation 

Central Vent 

Linear Vent 


Equation 



( 1 ) control volume 

V 

33 

il 

V = 2 bLdz 

(2) entrainment area 

A e 

A e = Inrdz 

A„ « 2 Ldz 

(3) mass: dry air 

~[p^ V( Pd\= pX<™)a c 

~-[p ll ur l <l> d ] = 2ap a ur 

- Upj ub 0ci] = a P« u 

az 

(4) mass: vapor 


^[ Pv ur\y 0 

~[p v ub<t> v \ = G 

(5) mass: solid 


= o 

~;[p,ub(p,} = 0 

(6) mass: bulk 

~[p B uV ] = p a (au)A e 

~\_PB url ] - 2 ap a ur 

^-{p B ub\ = ap a u 

(7) momentum 

~[p B u 2v ] = g(p a -p B )V 

~:[py'- 2 ] = g(p a -P B )r 2 

~[p B u 2 b] = g{p a - Pll )b 

(8) energy 

j-[p B uVC B 0] = 

^[p„ur 2 C B 0] = 

2ap u urC o e a - p o ur 2 g{^ d + <f > v ) 

~[p B ubC B 0] = 
ap a uC„9 a - p a ugb (f + <j > v ) 


542 

543 
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544 


Table 2. Notation’ 
Variable Definition 


A e 

b 

C a 

C B 

C v 

G 

g 

L 

n 

P 

r 

Rd 

Rv 

u 

V 

z 

a 

e 

Oa 

Pa 

PB 

Pd 


Pv 


Ps 


0s 

tv 


Cross Sectional Vent Area at the top of the conduit (m~) 

Entrainment Area (m 2 ) 

Linear Plume Half Width (m) 

Specific Heat of Ambient Atmosphere (Earth = 998 J K' 1 kg' 1 ; Venus = Mars = 
835 JK' 1 kg' 1 ) 

Bulk Specific Heat of Plume (J K' 1 kg' 1 ) 

Specific Heat of Water Vapor (= 2000 J K' 1 kg' 1 ) 

Specific Heat of Solid Particles (= 920 J K' 1 kg' 1 ) 

Gravitational Constant (Earth = 9.8 m s' 2 ; Venus = 8.4 m s' 2 ; Mars = 3.7 m s' 2 ) 
Fissure Length (m) 

Gas Mass Fraction (=p x <ffPn) 

Pressure (Pa) 

Cylindrical Plume Radius (m) 

Gas Constant for Ambient Atmosphere (Earth = 287 J K' 1 kg' 1 ; Venus = Mars = 
191 JK' 1 kg' 1 ) 

Gas Constant for Water Vapor (=461 J K' 1 kg' 1 ) 

Bulk Rise Velocity (m s- 1 ) 

Control Volume (m 3 ) 

Vertical Distance (m) 

Entrainment Constant (=0.09) 

Bulk Plume Temperature (K) 

Ambient Atmopshere Temperature (K) 

Ambient Atmospheric Density (kg m' 3 ) 

Bulk Plume Density (kg m' 3 ) 

Density of Entrained Ambient (kg m' 3 ) 

Density of Water Vapor (kg m' 3 ) 

Density of Solid Particles (= 1000 kg m' 3 ) 

Volume Fraction of Entrained Ambient 
Volume Fraction of Solid Particles 
Volume Fraction of Water Vapor 
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545 


o” subscript indicates a value at the vent, defined as the top of the conduit. 
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Table 3. Comparison of terrestrial vent sizes for circular and linear vents and a range of cross 
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Table 4. Comparison of volcanic vent sizes on 
cross sectional vent areas. 
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for circular and linear vents and a range of 


L 0 = 1000 

L 0 = 2000 

L 0 = 5000 

m 

m 

m 

0.628 

0.314 

0.126 

1.41 

0.707 

0.283 

3.93 

1.96 

0.785 

7.69 

3.85 

1.54 

15.7 

7.85 

3.14 

24.5 

12.3 

4.91 

35.3 

17.7 

7.07 

62.8 

31.4 

12.6 

79.5 

39.8 

15.9 

95.2 

49.1 

19.6 


70.7 


28.3 















































555 Table 5. Comparison of volcanic vent sizes on Mars for circular and linear vents and a range of 

556 cross sectional vent areas. 

557 


A 0 [lTf ] 

r» [m] 

b D [m] 

L 0 = 200 
m 

L 0 = 500 
m 

L 0 = 1 000 m 

L 0 = 2000 m 

L 0 = 5000 m 

849 

30 

2.12 

0.849 

0.424 

0.212 

0.085 

3150 

100 

7.89 

3.15 

1.57 

0.789 

0.315 

7150 

175 

17.9 

7.15 

3.57 

1.79 

0.715 

13700 

300 


13.7 

6.83 

3.42 

1.37 

25100 

500 


25.1 

12.5 

6.27 

2.51 

37800 

700 

■■ 

37.8 

18.9 

9.45 

3.78 

55700 

1000 

w 


27.8 

13.9 

5.57 

59900 

1200 

.v ,,': - 


30.0 

15.0 

5.99 
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